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ABSTRACT 

This paper constructs an analytic framework for calculating the assembly of 
galactic disks from the collapse of gas within dark matter halos, with the goal 
of determining the resulting surface density profiles. In this formulation of the 
problem, gas parcels (baryons) fall through the potentials of dark matter halos on 
nearly ballistic, zero energy orbits and collect in a rotating disk. The dark matter 
halos have a nearly universal form, as determined previously through numerical 
simulations. The simplest scenario starts with a gaseous sphere in slow uniform 
rotation, follows its subsequent collapse, and determines the surface density of 
the disk. This calculation is carried out for pre-collapse mass distributions of 
the form M(r) ~ with I < p < 3 and for polytropic spheres in hydrostatic 
equilibrium with the halo potential. The resulting disk surface density profiles 
have nearly power-law forms a ~ w^'^ (where w is the radial coordinate in the 
disk), with well-defined edges determined by the centrifugal barrier Rc - the 
radius to which gas with the highest specific angular momentum falls during the 
collapse. This idealized scenario is generalized to include non-spherical starting 
states, alternate rotation profiles, and multiple accretion events (e.g., due to gas 
being added to the halo via merger events). This latter complication is explored 
in some detail and considers a log-normal distribution for the angular momenta 
of the pre-collapse states of the individual components. In the limit where this 
distribution is wide, the composite surface density approaches a universal form 
(Tt ~ ti7~^, independent of the shape of the constituent profiles. When the angu- 
lar momentum distribution has an intermediate width, however, the composite 
surface density attains a nearly exponential form, roughly consistent with profiles 
of observed galaxies. 

Subject headings: galaxies: formation, halos, kinematics and dynamics — meth- 
ods: analytical 
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1. INTRODUCTION 

Although the formation of dark matter halos is rapidly becoming understood, the for- 
mation of gaseous galactic disks remains more elusive. Disk formation is complicated and 
involves a number of physical processes, including the large scale collapse of the dark mat- 
ter, merging of separate halos, gas cooling, and the infall of gas through the gravitational 
potential of the dark matter. Numerical simulations now include all of these processes and 
are beginning to produce reahstic gaseous disks (e.g., Governato et al. 2004; Robertson et 
al. 2004; Keres et al. 2005). This paper adopts a complementary approach and seeks to 
study the infall process — separately and analytically — where the starting conditions are 
informed by numerical simulations that have been previously carried out. This paper thus 
develops an analytic framework from which to study the infall process in galactic disk forma- 
tion. This formalism applies to a wide range of infall scenarios and can be readily generalized 
to include additional physical processes. The main focus of this paper is to investigate the 
range of surface density profiles produced by gaseous collapse within the potential of a dark 
matter halo. 

Although the aforementioned simulations demonstrate that the current paradigm of 
structure formation can produce galactic disks with reasonable properties, an analytic de- 
scription of disk formation is useful in several respects. Analytic solutions provide us with a 
direct assessment of how the final disk surface density properties depend on starting condi- 
tions and other parameters of the problem. Along the way, we obtain an analytic description 
of several ancillary quantities, e.g., the velocity and density fields of the gaseous material 
as it falls toward the disk, as well as the hydrostatic equilibrium profiles of the pre-collapse 
states. Once found, the resulting analytic solutions, both for the ancillary quantities and 
the resulting disk surface density profiles, can be used in a variety of other contexts. 

In the standard picture of galactic disk formation (e.g., Rees & Ostriker 1977; Fall & 
Efstathiou 1980; White & Frenk 1991), the gas within a forming dark matter halo approaches 
a density profile that is roughly parallel to that of the dark matter. This assumption is jus- 
tified by complementary simulations of gas dynamics (from Evrard et al. 1994 to Kaufmann 
et al. 2006). In order for the gas to fall inward toward the nascent galactic disk, the gas 
must cool. The gas cooling time depends primarily on density, which is a decreasing function 
of galactic radius, so the cooling time increases outward (note that metallicity, temperature, 
etc., also influence the coohng time). There are thus two regimes of interest. In one hmit, 
the entire virialized region has a cooling time shorter than the age of the universe (at the 
formation time of the halo). In this case, the gas cools rapidly and never approaches hydro- 
static equilibrium. The gas supply is then limited by the rate at which the forming galactic 
structure accretes gas from larger radii. In the other limit, the cooling time is longer than 
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the age of the universe for much of the viriahzed region, and infall is hmited by the rate of 
gas coohng (White & Rees 1978). The flow of gas into the central regions is then analogous 
to cooling flows (e.g., Fabian et al. 1984). In both regimes, the gas falls inward through the 
gravitational potential provided (primarily) by the dark matter halo. 

In this paper, we study the infall of baryonic gas moving in the potential of a dark matter 
halo. Since dark matter halos approach a known nearly universal profile (e.g., Navarro et al. 
1997; hereafter NFW), we consider the potential to have a given form and follow the orbits 
of gas parcels as they fall through the dark matter and build up a galactic disk structure. 
Specifically, wc consider the halo to have a fixed profile, independent of the gas, and to 
take the form of a Hernquist potential (Hernquist 1990; eq. [1]). Notice that this form is 
essentially the same as the more commonly used NFW profile, except that it has a steeper 
density profile at large radii and hence converges to a finite mass (see also Busha et al. 2005). 
As outhned above, the gas must cool before it falls through the halo potential toward the 
central regions. For the first part of this work, we describe this process as a single, overall 
collapse, even though the gas cools at different rates. This complication will not affect the 
resulting surface density profiles of the disk as long as the central portions fall in first (as 
expected, since they tend to be denser and cool faster) so that infalling shells do not cross. 
Notice also that we are implicitly assuming that the gas mass is small compared that of the 
dark matter. After a sufficiently large quantity of baryonic mass has accumulated near the 
galactic center, however, the back reaction of the gas on the dark matter potential should 
be included. 

The goal of this calculation is thus relatively modest: We want to develop an analytic 
transformation between starting gas distributions and the "final" surface density profiles of 
the disk. In this context, the starting gas distributions refiect the initial profiles of specific 
angular momentum, which depend on gas cooling and other effects that determine the start- 
ing state. The final surface density profiles (as considered here) are determined by where gas 
falls onto the galactic disk with no additional evolution (e.g., due to gravitational torques, 
which act on longer time scales). We discuss the effects of additional evolution, e.g., angular 
momentum redistribution, in the final section. 

This work expands upon previous work concerning galactic disk formation (sec, e.g., 
Gunn 1982, Kauffmann 1996, Dalcanton et al. 1997, Firmani & Avila-Reese 2000, Bullock 
et al. 2001, van den Bosch 2001, Thacker & Couchman 2001). These previous studies ex- 
plore disk formation within dark matter halos with increasing levels of complexity, primarily 
through numerical simulations. The generic finding is that the resulting surface density 
profiles are more concentrated in the center than an exponential law, and typically display 
sharper outer boundaries. This result depends on the starting (pre-coUapse) distribution 
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of angular momentum of the gas. Dalcanton et al. (1997) start with a uniformly rotating 
sphere, Firmani & Avila-Reese (2000) assume that each infalling mass shell has a constant 
spin parameter A, and Bullock et al. (2001) determine the initial angular momentum distri- 
bution directly from N-body simulations. In these studies, the disk surface density profiles 
are the result of angular momentum conservation, although the initial distributions of an- 
gular momentum vary. One contribution of this present study is to provide a systematic 
exploration of possible starting states and their effect on the corresponding surface density 
profiles. In particular, we generalize the starting conditions to include a wide variety of 
initial density profiles, rotation profiles, and initial geometries. We also explicitly determine 
the velocity fields and density fields of the infalling material. More significantly, this analytic 
framework is sufficiently robust to account for multiple accretion events (see also Firmani 
& Avila-Reese 2000). For comparison, previous numerical work (Robertson et al. 2004; 
Keres et al. 2005) includes more physical processes (e.g., explicit calculation of gas cool- 
ing and angular momentum transport during collapse), but does not isolate the dynamical 
mechanism(s) by which galactic disks are assembled. 

This paper is organized as follows. We first extract the relevant orbit solutions that 
describe gas parcels falling through the gravitational potential of a dark matter halo (§2); 
this calculation determines the velocity and density fields of the infalling material, as well 
as the locations where given parcels fall onto the disk plane. In §3 we apply these results 
to the collapse of rotating gaseous spheres and find the resulting surface density profiles as 
a function of the starting density and angular momentum distributions. These results are 
generalized (in §4) to include more complex starting geometries, including holes, removal of 
polar caps, and filamentary structures. Motivated by the fact that cosmic structure is built 
up through a hierarchical process, §5 considers the effects of multiple accretion events and 
finds the corresponding composite surface density profiles. The paper concludes in §6 with 
a summary of our results and a discussion of how exponential disks can be produced. 



2. DYNAMIC INFALL IN A DARK MATTER HALO POTENTIAL 

For a given gravitational potential, we must find orbital solutions for gaseous material 
falling towards the galactic center. Here we assume that the halo has the form of a Hernquist 
profile so that the potential, density distribution, and mass profile are given by 

where ^ — r/vs and where Vg is the scale length of the structure. We expect ~ 50 — 80 
kpc for a typical galactic halo. The other constants are related to each other through the 
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expressions = ^nGpoVg and Moo — 27rrfpo- Notice that is defined to be positive (so 
that the signs in the equations of motion must be chosen accordingly). Notice also that 
the halo parameters (\l/o, r^, Mqo) will evolve with time and will, in general, be different for 
different accretion events. As a starting approximation, however, these quantities are taken 
to be constant during a given accretion event. 

We must specify the initial conditions for the infalling baryonic gas. For any orbit, the 
initial condition can be specified by the energy and angular momentum. As first approxi- 
mation, we consider the gas parcels to enter the central region of the galaxy on zero energy 
orbits. This assumption could be relaxed to include small but nonzero energies, which can 
be written in dimcnsionlcss form 



In addition, the gas parcels are assumed to have a given specific angular momentum j — j^o 
at their starting radius. The angular velocity vector of the starting (spherical) state points 
along the z axis, but the angular momentum vector of an individual gas parcel depends 
on its polar angle (see below). In order for the gas to fall relatively close to the galactic 
center - onto the forming galactic disk - this angular momentum must be sufficiently small. 
Following previous work (Adams & Bloch 2005; hereafter AB05), the orbits are described 
by the angular momentum parameter 



The condition that the orbits fall in the central region of the halo thus takes the form <S 1. 
In the limit of interest, -C 1 and e — > 0, the orbit solutions take the approximate form 



where is the turning angle in the plane of the orbit (AB05) . Since the potential is spherically 
symmetric, angular momentum is conserved and the motion is confined to a plane described 
by the coordinates (r, 0) ; the radius r is the same in both the plane and the original spherical 
coordinates. The angular coordinate in the orbital plane is not the same as the polar angle 
6 in the original spherical coordinates, but the two are simply related. Suppose that the 
gas parcel has polar angle 6q at the start of its orbit, i.e., we let 6q denote the angle of the 
asymptotically radial streamline. The coordinate position of the gas parcel is subsequently 
described by either the angle in the orbital plane, or by the polar angle 9 in the original 
spherical coordinates (see Figure 1). The two angles obey a geometrical relation (e.g,. Ulrich 
1976) that depends on the initial polar angle, i.e.. 



e= \E\/^o < 1. 



(2) 



(3) 



^ sin = ^/rj — constant , 



(4) 



COS0 = 



COS^ _ fjL 



(5) 



cos^o A*o ' 
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where the second equality defines fi = cos 9 and no = cosOq. Combining the above relations, 
we can write the orbit equation (4) in terms of the variables /i and /Iq to obtain 

For zero energy orbits, conservation of energy enforces a constraint on the velocity fields 
of the form 

-r+vl+vl^^^. (7) 

Given the orbital solution, we can find the velocity fields 

(8) 

Since ^ and /iq are related through the orbit equation (6), the velocity field is completely 
determined for any position (r, 6) . The density distribution of the infalling material can 
be obtained by applying conservation of mass along a stream-tube (Ulrich 1976; Chevalier 
1983), i.e., 

p{r, 9) Vr sin 9 dO dip = — — sin 9o d9o dipo, (9) 

where M is the total rate of mass flow through a spherical surface (note that M is defined in 
the intermediate asymptotic regime inside the starting radii and outside the galactic disk). 
The minus sign arises because the radial velocity Vr is negative for infalling trajectories. 
Combining the above equations allows us to write the density profile of the incoming material 
in the form 

M duo 

The properties of the collapsing structure determine the orbit equation (6), which in turn 
determines the form of djiQ/dji (see also Cassen & Moosman 1981; Terebey et al. 1984). 

With the density specified, we can determine the rate at which the forming galactic 
disk gains mass, da/dt — pve, where the right hand side of the equation is evaluated at the 
disk plane 9 — 7r/2 or /i — 0. Notice that we are assuming that the newly arriving gaseous 
material dissipates the energy contained within its f and 9 velocity components, so that it 
joins the disk. Such dissipation can take place, e.g., by shocks with the pre-existing disk 
material, or, in the case of sufficiently symmetric flow, by meeting the corresponding gas 
parcel arriving at the disk plane from the opposite direction. With these assumptions, the 
disk grows according to 

da - — — (11) 
47rr^ Vr djJL 



Fig. 1. — Schematic diagram showing the geometry of the orbits. The heavy sohd curve 
depicts the orbit trajectory, which begins at large radius tq ^ r and continues inward until 
it reaches the disk plane. The orbit proceeds within a fixed plane that lies at an angle 9q, 
measured with respect to the z axis {tq-z = cos^^o)- At a given point along the orbit, marked 
by position vector r = rf, the angle within the orbital plane is (r ■ fo = cos0) and the 
position angle in the original coordinate system is ^ (r • i = cos 9) . The thick fine marked 
by square symbols shows the intersection of the orbit plane with the disk plane. 
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Notice that the above equation only accounts for the surface density distribution pro- 
duced by the arrival of the gas. After its arrival, the gas will dissipate further energy, transfer 
angular momentum, and the disk will evolve. In particular, the newly arriving material does 
not know in advance that it will become part of a centrifugally supported disk structure. As 
a result, it will not necessarily have the proper azimuthal velocity appropriate for the local 
circular speed of the disk (e.g., Cassen & Moosman 1981). As it adjusts to disk conditions, 
the gas will dissipate additional energy and move deeper into the gravitational potential 
well of the galaxy. On longer time scales, the disk will experience gravitational torques, 
which transfer angular momentum and act to spread out the disk (e.g., Shu 1990; Binney & 
Tremaine 1987). These effects will be discussed in greater detail below (see §6.2). 



3. BASIC COLLAPSE SOLUTIONS 

In the standard picture of galactic disk formation, baryonic matter initially traces the 
density distribution of the dark matter halo and attains a similar distribution of angular 
momentum. After cooling, the gas falls inward on nearly ballistic (pressure-free) trajectories 
and conserves its angular momentum. The objective of this work is to understand this 
infall process in greater detail, starting form simple models and then generalizing to more 
realistic models. In this section we consider the surface density profiles produced through 
the collapse of a single, well-defined gaseous configuration embedded within the potential of 
the dark matter halo. We consider more general collapse scenarios, including non-spherical 
effects and multiple accretion events, in the following sections. 

To start, we must thus specify the starting mass profile of the (baryonic) gas. Here we 
assume that the mass profile M(r) is spherically symmetric and that the inner regions fall 
before the outer regions so that infalling shells do not cross. We also assume that the initial 
structure is slowly rotating with a well-defined rotation profile fl{r). Note that when a given 
mass M has fallen to the disk, the inverse relation r(M) determines the initial location - and 
hence the angular momentum - of the newly arriving material. The mass profile M(r) thus 
helps set the angular momentum profile of the starting state for a given il(r), and thereby 
influences the resulting surface density proflle. 

If the gas has an initially spherically symmetric density distribution and is rotating with 
angular velocity ft, as described above, the specific angular momentum for a given starting 
position (^co, ^o) is given by 

ioo^r^e^^^sin^o, (12) 
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which implies that the angular momentum parameter takes the form 



^=^d(l-/^o)- (13) 



The orbit equation (6) thus becomes 



1-^^ = 1(1-/^0), (14) 

where we have defined 

^0 - "-^C - , (15) 

where the second equality defines the dimensionless rotation parameter cu. 

With the above considerations, the density field is thus completely determined in ana- 
lytic but implicit form 

M l-i^l 

Notice that the density must be expressed in terms of spherical coordinates r and 6 = cos~^ fi, 
but the right hand side of the equation is written in terms of /^o = cos^o- The two variables 
arc related through the orbit equation (14). 



P(e,^) = 7-r^^- (16) 



For the case of spherical collapse, we can now evaluate d/jio/dfi, the velocity fields, the 
density profile, and hence the rate at which the disk gains mass from the infall. We thus 
obtain 

da = e = 7r/2)vedt = - ^7%]"'^' , (17) 

where we have multiplied the right hand side by a factor of 2 to account for gas fiowing 
onto the disk from both sides (top and bottom). In this model, tjq oc so the the starting 
mass profiles [basically, how the enclosed mass M(r) depends on r] determine how much 
mass initially has a given specific angular momentum and hence the radial location on 
the disk where it will fall. Notice that the derivation of equation (17) uses an extreme 
approximation in which the turning point of the orbit occurs at the point = 0, i.e., where 
the orbit intersects the disk. As a result, both the radial velocity Vr and the differential 
d/io/d/i contain integrable singularities, which must be handled properly. As expected, these 
singularities cancel out and reveal a finite mass fiow onto the disk. 

For each given mass and angular momentum profile, the above formalism allows us to 
calculate the corresponding surface density of the resulting disk (see §3.1 - §3.5). We begin 
with some general definitions. For any given mass profile, one can define a mass scale Ms to 
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be the mass in the gaseous sphere contained within the scale radius of the dark matter 
halo, i.e., 

= M(r,) . (18) 

This mass scale is a convenient reference mass. Notice that Mg is the baryonic mass and 
is thus a small fraction (~ 1/6, Spergel et al. 2003) of the total mass contained within the 
radius r^. We can also define a rotational r^^ within the collapsing gaseous sphere, i.e., 

r^ = rsu; = rlQ/{2^of'. (19) 

As shown below, the centrifugal barrier, which sets the location of the outer disk edge, is 
given by a dimensionless factor times the scale r^^. As gaseous material falls inward, it has a 
given radius r in spherical coordinates, and then becomes part of a disk which is generally 
described in terms of cylindrical coordinates. The radius r = r^^ in spherical coordinates 
thus becomes the radius within the disk in cylindrical coordinates and is denoted as 

w = rs^. (20) 

With these definitions, we can find the surface density profiles for a given starting state. 



3.1. Specific spherical starting states with constant rotation rate 

We can consider a range of specific models for the starting mass profile, including a 
self-gravitating isothermal sphere 

, , , , 2a2 , ^ ^ GM , , 

pW = ^^, M{r) = —r, and r(M) = ^ , (21) 

where a is the isothermal sound speed. We also consider a uniform density sphere 

p(r)=po, M(r) = ^por^ and K^) = (|^) , (22) 

and an intermediate density profile motivated by the Hernquist (and NFW) profile, where 
the (pre-coUapse) density profile takes the form 

p(r) = Po(ro/r) , , M(r) = 2wor' , and r(M) = . (23) 

This profile would occur, e.g., if the baryonic gas traces the density distribution of the dark 
matter halo in the inner regime. For this profile, without loss of generality, we can take Tq 
= Tg and define the constant po accordingly. 
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Starting with the case of the isothermal sphere (eq. [21]), we find 

This expression can be written in terms of the dimensionless variable x = (M/Ms)(a;/4)^/^, 
and hence the integral expression for the surface density becomes 

The range of integration starts at a; = 1, rather than x = 0, because all of the material 
with X < 1 falls inside of the radial coordinate w in the disk. For a given physical radius w 
within the disk, the innermost portion of the sphere has too little angular momentum (even 
at the equator where angular momentum is maximum) to fall as far out in the disk as the 
radial location w. The upper limit of integration is set by the total mass Mt of the accretion 
event, i.e., Xf = [Mt / Ms){r^ / wY^"^ . In this scenario, the disk has a well-defined outer edge 
determined by the material with the highest specific angular momentum in the initial state. 
We denote this radius as the centrifugal barrier Rc, which takes the form Rq — r^{MT/MsY, 
so that Xf — {Rc/voY^'^. The surface density can thus be written in the form 

( \ pc^^ri -41-1/2 /ofi^ 

Note that when w ^ Rc-, far from the outer disk edge, the surface density attains a nearly 
power-law form (here a ~ -07"^/^). Near the edge, where zu — >■ Rc-, the upper limit of the 
integral Xf 1 and the surface density a — > 0. 

A similar treatment can be used for the other mass profiles. For the case of a uniform 
density sphere (eq. [22]), we define x = {M/Ms){u!/^Y^'^, and find the surface density 

where Rc = r^^^Mx / Mg)"^^^ and hence Xf = {Rc/'c^Y^'^- For the intermediate case (eq. [23]), 
we define x = {M / Ms){lo / and find the surface density 

aivj) = / —\l—x ' = — cos ivj Rc) , (28) 

where Mt is the total mass that has fallen to the disk and Rc — {Mt/ Ms)ruj. 

To illustrate the surface density profiles that result from a single spherical collapse. 
Figure 2 shows the results for three mass profiles M(r) ~ r^, where p = 1, 2, and 3. All 
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Fig. 2. — Surface density profiles resulting from spherical collapse for three values of the 
index p of the starting mass profile M(r) (x r^: p — 1 (dotted curve), p — 2 (sohd curve), 
and p — 3 (dot-dashed curve). The rotation rate of the starting states is taken to be a 
constant. All three profiles have the same total mass and the same value of the centrifugal 
barrier (scaled here so that Rc — 24 kpc), which determines the outer disk edge. The dashed 
curve shows an exponential profile for comparison, where the exponential scale length is i?x 
= 5 kpc and the total disk mass is the same as the other profiles. The scale on the vertical 
axis is arbitrary. 
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three profiles have the same total disk mass and the same centrifugal radius, which sets the 
location of the outer disk edge. These profiles are scaled so that Rc = 24 kpc, but solutions 
can be found for any disk mass and size Rq- To fix ideas, consider halo/disk parameters 
roughly comparable to those of the Galaxy. If we use a Hernquist profile to describe the 
Galactic potential, the scale length ^ 65 kpc (see AB05 for further detail) and the mass 
profile for the gas corresponds to p = 2 if we invoke an outer boundary at rr ~ A galactic 
disk with Rc — 24 kpc would require an initial rotation rate Q ~ 10~^^ sec~^ and would 
result in a disk mass ~ IO^^Mq. This value of Q is comparable to (but somewhat larger 
than) that expected for a spin parameter A 0.05 (Bullock et al. 2001) for the same halo. 
This disk mass implies a total galactic mass of about 2.4 x 10^^ Mq (since only 1/4 of the gas 
mass collapses and the ratio of baryons to total matter is about 1:6).^ Figure 2 also shows 
an exponential surface density profile for comparison, where the scale length R^ = 5 kpc 
and the exponential disk has the same total mass. Compared to the exponential disk, the 
surface density profiles resulting from these spherical collapse solutions increase more steeply 
toward the center and truncate sharply at the edge {Rc)- Nonetheless, over the intermediate 
range of radii shown here, the surface density profiles are roughly similar. Notice also that 
these profiles are similar to those obtained in previous work based on conservation of angular 
momentum only (see Fig. 1 of Dalcanton et al. 1997) and by detailed SPH simulations of 
disk formation including many additional physical processes (compare with Figs. 9, 10, and 
16 of Kaufmann et al. 2006, and with Figs. 10 and 11 of van den Bosch 2001). 

Although we generalize this basic calculation of the surface density in subsequent sec- 
tions, the general features shown in Figure 2 are robust, namely a power-law profile with a 
relatively sharp truncation at the outer disk edge. This outer edge is set by conservation 
of angular momentum. For a given starting state, the gas has a maximum specific angular 
momentum jmax- Material with angular momentum j = j^ax will fall to a well-defined radius 
(that of the disk edge) and material with j < j^ax will fall to smaller radii. The surface 
density profiles tend to diverge in the center — > 0. The centers of actual galaxies are 
complicated by the presence of bulges (Binney & Merrifield 1998) and central black holes 
(Gebhardt et al. 2000). Further, this behavior can be controlled by including a core radius 
in the initial density profile (§4). 

In order to understand the infall picture in greater detail, it is useful to see where gas 
parcels that start at various initial locations end up within the disk. Figure 3 illustrates 
the general behavior (using the p — 2 mass profile). The surface density at a given radial 



For completeness we note that the baryon to dark matter ratio in galaxies is lower than that of the 
universe as a whole by a factor of 3 5. This difference can be explained by internal galaxy processes, e.g., 
long gas cooling times coupled with strong feedback. 



-14- 



location within the disk, denoted by the coordinate ro, is fed by a range of initial radii r 
within the starting (spherical) state. As the initial radius increases, the initial polar angle 
9 decreases along a well-defined locus of points, as shown in the left panel of Figure 3. For 
a given radial location w in the disk, there exists a minimum polar angle that contributes 
to the surface density; for starting polar angles constrained by the relation sin^ < w/Rc-, 
all the of material falls to smaller radii w' < w. One can also ask where the material at a 
given starting radius ends up. This result is shown by the right panel of Figure 3, where 
gas parcels along the outer boundary of the collapsing region are mapped to their landing 
points on the galactic disk. The pole {9 — 0) is always mapped to the galactic center and the 
material along the equator falls to the outer disk edge (in this geometry, the material along 
the equator has the highest specific angular momentum and defines the outer disk edge). 



3.2. General spherical starting states with constant rotation rate 

The treatment developed here can be generalized to include any power-law mass profile, 

i.e., 

M(r) oc and r(M) oc M^'^ . (29) 

We define the dimensionless variable x = {M/Ms){uj/^y/'^, and the surface density profile 
can be written in the form 



where 

i?c = r,(Mr/M,)2/^ (31) 

so that Xf = {Rc/tu)p/^. 

As a consistency check, we can show that the total mass locked up in the surface 
density profile is equal to the original mass Mt that falls onto the disk. To show this, we 
integrate the general form of the surface density (eq. [30]) over the surface of the disk to 
obtain 

Md = 27rwdw ,^ r ^[1 - x-^'P]-^''' . (32) 

The outer edge of the disk is given by Rc — r^{MT/Msy/P and the upper limit of integration 
in the first integral can be written as Xf — {Rc / vjY^^ . If we change the integration variable 
from zu to z = w/Rc, the integral becomes 



(V^)^/^ dx_ 



Mn-Mrl-^J - x-^'r'" ■ (33) 
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Fig. 3. — Diagram showing the relationship between starting locations for the orbits and the 
final locations. The left panel shows the starting positions that lead to gas parcels falling to 
a given radial location within the disk; the dotted curve shows the locus of starting points 
for which the orbits land at the particular location -07 = 10 kpc. The solid curve shows the 
outer boundary of the collapsing region (taken here to be 40 kpc). The right panel shows 
the final radial positions in the disk for starting points on a given spherical surface, taken 
here to be the outer boundary at 40 kpc. In both panels, the dashed fines are schematic and 
connect the starting and ending locations. 
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Next, we change the order of integration and evaluate: 

t J^i---r-r''';^^^f ^-ii-.-r- (3.) 

This argument shows that mass is conserved for any mass profile, as characterized by the 
index p, i.e., the resulting surface density profile has the same total mass as the starting 
state. 

The general result of this analysis is that a single spherical collapse naturally produces 
a power-law surface density profile with a sharp truncation at the outer edge. We can be 
more precise about the sharpness of the edge: For any value of p, we obtain a integral / 
(see eq. [30]) that is a function of Xf and hence a function of the radial disk coordinate zu. 
The variation of the surface density profile can be characterized by the index Q defined by 
Q = {w / I){dl / dw). It is straightforward to show that Q ^ oo in the limit that vo — > Rc 
for all indices p > 0. 



3.3. Starting states with constant rotation speed 

Many numerical simulations of structure formation suggest that the angular momentum 
distribution of the gas, before collapse, takes the form j(r) oc r (Kaufmann et al. 2006; 
Bullock et al. 2001). This angular momentum profile is equivalent to a rotation rate with 
radial dependence oc r^^, rather than the simplest case of ^2 = constant considered above. 
We can generalize the previous treatment as follows. The rotation rate of the starting gas 
can be written in the form 

n{r) = Qsj = , (35) 

where we have scaled the angular velocity to its value at r = r^. For this rotation profile, 
the angular momentum parameter r]Q (see eq. [15]) takes the form 

* = g?i = -'fi where . ^ . (36) 

For a given mass profile of the form M oc r^, the integral expression that determines the 
surface density profile takes the form 

where the integration variable is given by a; = {M/Ms){ri^/zuy. For these models, the 
centrifugal radius, which sets the location of the outer disk edge, is given by 

Rc^r^Mr/M^y/P, (38) 
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and hence Xf — {Rc/vjy. 

For this choice of initial rotation profile, the surface density integrals can be evaluated 
in terms of elementary functions for starting mass profiles with index p = 1, 2, and 3. For 
the starting density profile of an isothermal sphere, p = 1, we obtain 

= 2^^^ cos-\w/Rc) , (39) 

which is the same as that obtained earlier for the case of the intermediate p = 2 mass profile 
and a constant rotation rate. With this rotation profile VL ~ r^^ and p = 2, the surface 
density takes the form 

Mr 

a{zu) = cos\i-\Rc/vj) . (40) 
For a uniform density starting state, p = 3, the surface density becomes 



-(-) = |^(1--Vi?^)^/^ (41) 

Note that equations (39 - 41) describe how the surface density profiles (for constant rotation 
speed) depend on the underlying parameters (where the definitions of eqs. [36] and [38] 
complete the specification). 

Figure 4 shows the surface density profiles for the three mass profiles p = 1, 2, and 3, 
for the case of constant starting rotation speed [Vt ~ r~^). These profiles arc less steep than 
those with constant rotation rate (fi = constant] see Fig. 2). In this case, thep = 1 starting 
state leads to a power-law surface density profile a ~ zu~^ and the p = 3 starting states 
leads to (T ~ (where both power-laws are truncated at Rc by an edge function). The 
intermediate mass profile with p = 2 produces a surface density profile of the approximate 
form a ~ log(i?c/'Ci7) away from the disk edge (i.e., for zu <^ Rc)- For ~ a vahic of 
the constant rotation speed Vg — Tgfls ~ 180 — 200 km/s is required to produce Rc — 24 kpc 
as shown here. 



3.4. General spherical starting states 

This subsection considers a general spherical starting state in which the mass profile has 
the form M (x (with p in the range 1 < p < 3, but otherwise arbitrary) and the rotation 
profile has the form fl oc (with s in the range < s < 1). The centrifugal radius takes 
the form 

/ Mt \ ^/"^ p 
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Fig. 4. — Surface density profiles resulting from spherical collapse with starting rotation rate 
Q oc r"^. The resulting surface density is shown for three values of the index p of the starting 
mass profile: p = 1 (dotted curve), p = 2 (solid curve), and p = 3 (dot-dashed curve). All 
three profiles have the same total mass and the same value of the centrifugal barrier (Re 
— 24 kpc), which determines the outer disk edge. The dashed curve shows an exponential 
profile for comparison, where the exponential scale length is i?x = 5 kpc and the total disk 
mass is the same as the other profiles. The scale on the vertical axis is arbitrary. 
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The index a plays an important role in determining the functional form of the surface density, 
which is given by 

a{w) = r ^[1 - . (43) 

The effective power-law index q of the surface density profile (far from the disk edge at Rc) 
is given by 

g = 2 - a = 2 - ^ . (44) 

This interpretation (that q is the power-law index of the profile) breaks down when q ^ 0, 
i.e., the portion of parameter space defined by g < results in qualitatively different surface 
density profiles (ones that cannot be described in terms of power- laws). In general, these 
surface density profiles are shallower than power-laws. 



3.5. Hydrostatic initial states 

This subsection considers the scenario in which gas enters the central region of a dark 
matter halo and reaches hydrostatic equilibrium before cooling and falling inward to the 
disk. In this approximation, the gravitational potential is provided by the dark matter halo 
alone. The resulting (pre-coUapse) density profile depends on the equation of state for the 
gas. We first consider the idealized case of an isothermal gas with a uniform rotation rate; 
we then consider the more realistic case of a general polytropic equation of state with a 
constant velocity (O oc 1/r) initial condition. 

In the limit of an isothermal equation of state, the starting hydrostatic density profile 
takes the form 



■(^o/a2)i 



P(r) = Poo exp [———] = Po exp -[^) — — , (45) 



1 + e ' """""L VaV 1 + e 
where a is the isothermal sound speed of the gas. After coohng, the gas falls through the 
potential of the dark matter halo and collects into a disk structure. In this case, the integral 
that determines the disk surface density is more naturally written as a radial integral, rather 
than a mass integral. We assume that the starting (pre-coUapse) density profile follows the 
form of equation (45) out to an outer boundary radius tt- We then work in terms of the 
variable defined by y = (r/rs)(r^;/ti7)^/^, and define the centrifugal radius via 

Rc^rUrr/rsf. (46) 

Note that we have assumed a uniform rotation rate in deriving this form for Rc- The integral 
that determines the disk surface density then takes the form 

2rf,po f^f dy r /*o\ Ay i 
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where A = {w /RcY''^{rT/rs). 

Figure 5 shows the surface density distributions resulting from this type of collapse, 
where the integral in equation (47) has been evaluated numerically. For the cases shown 
here, the initial radial extent of the pre-coUapse gas was taken to be tt — Ts, the rotation 
rate was chosen so that the centrifugal radius Rc — 24 kpc (i.e., Q ~ 10"^^ rad~^ for the 
Galactic parameters used previously), and all profiles have the same total mass. The four 
curves shown in Figure 5 correspond to varying values of the ratio "^o/a"^ = 2, 4, 8, and 
16 (from top to bottom in the figure). For moderate values of the ratio "^q/o^, the surface 
density distributions have the same general form as those obtained previously (compare with 
Figs. 2 and 4). 

This approach can be generalized to include any polytropic equation of state for the 
initially hydrostatic gas and any rotation profile fl{r). We start by writing the equation of 
state in the form 

where n is the polytropic index (e.g., Chandrasekhar 1939). The density distribution for 
hydrostatic equilibrium, in the limit where the dark matter halo dominates the gravitational 
potential, then takes the form 

P-T-^^ (49) 
(1 + 0" 

where the central density is specified through the relation 

po = [*on/(n + l)«]". (50) 

In the limit n ^ oo, the equation of state approaches an isothermal form; in the limit 
n— >3/2,7— >5/3, and the equation of state becomes adiabatic (for a monatomic gas). The 
polytropic index is thus confined to the range 3/2 < n < oo. We can see immediately that 
the value n = 3 (7 = 4/3) corresponds to a critical value: If n > 3, then the mass integral 
converges, whereas the mass integral diverges for n < 3. 

Following the same formulation used above for the isothermal case, the disk surface 
density resulting from the collapse of a polytropic sphere with rotation profile Q oc 1/r takes 
the form 

, , 2poTlw f" ydy 1 

= - 1)'" (1 + A,)» ■ ^''^ 

where i/f — {rT/rs){r^^/w) and A = [w |r^^. As a consistency check, we can integrate this 
expression over the surface of the disk to find the total disk mass, and thereby obtain the 
expected result 

Md = 27rwdwa{w) - inpor^ Jl+W " ^^''^^ ' ^^^^ 
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Fig. 5. — Surface density profiles resulting from spherical collapse starting from quasi-static 
equilibrium of an isothermal gas in a dark matter halo. The four curves show the surface 
density profiles for varying values of the ratio '^o/a^ = 2, 4, 8, and 16, where a is the isother- 
mal sound speed of the gas (before it cools and collapses). The surface density distribution 
becomes steeper as the ratio "^o/a^ increases. Cold gas (with sound speed a small compared 
to the depth of the gravitational potential well) results in a steep surface density profile be- 
cause mass is concentrated at small radii in the hydrostatic (pre-coUapse) state. The dashed 
curve shows an exponential profile for reference, where the scale length Ry^ — 5 kpc. The 
scale on the vertical axis is arbitrary. 
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Fig. 6. — Surface density profiles resulting from spherical collapse starting from quasi-static 
equilibrium in a dark matter halo where the gas has a polytropic equation of state and 
the starting rotation profile Q(r) oc 1/r. The profiles shown in the figure correspond to 
polytropic indices n — 2 - 16 (from top to bottom). All of the profiles are normahzed to 
have the same mass. As the polytropic index n increases, the equation of state becomes 
closer to isothermal, more gas is concentrated at small radii in the pre-coUapse state, and 
the disk surface density becomes steeper. The dashed curve shows an exponential profile 
with scale length = 5 kpc. The slope of the surface density produced from the collapse 
of the n = 8 polytrope is close to this reference profile. The scale on the vertical axis is 
arbitrary. 
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For any polytropic index n, the resulting surface density distributions display the same 
general form - with radial dependence somewhere between a truncated power-law and an 
exponential. In these disks, the outer radius is determined by the centrifugal barrier, which 
in turn is determined by the material with the highest specific angular momentum, and that 
material comes from the equatorial region at the outer boundary tt of the collapsing region. 
A collection of surface density profiles of this type is shown in Figure 6 for polytropic indices 
n — 2 - IQ. For the profiles shown here, the central density po is assumed to be constant and 
the profiles are normalized to have constant total mass. For these parameter choices, the 
surface density profiles become steeper as the polytropic index n increases. The resulting 
curves are close to straight lines on the plot of logcr versus w, where the effective scale 
radius of the exponential decreases with increasing polytropic index n. Figure 6 also shows 
an exponential profile with scale length = 5 kpc (dashed curve). The surface density 
profile produced by the collapse of the n = 8 polytrope has nearly the same slope. For 
comparison, models of galaxy groups and clusters have been used to estimate the equation 
of state in astrophysical systems (Ascasibar et al. 2003) and find a fit of the form P oc p^'^^, 
which corresponds to a polytropic index n ~ 5.6. 

For completeness we note that for all polytropic indices n > 3, the total mass of the 
pre-coUapse gaseous sphere is finite and disk solutions exist for the limit y/ — > oo. In this 
case, the resulting surface density profiles approach power-law forms such that a ~ zu'^"'^^ 
at large radii. However, the convergence to this power-law form is painfully slow. 

3.6. Section summciry 

Actual galactic disks differ from the idealized disks found in this section in several ways. 
As outlined below, the single, spherical collapse picture must be generalized to allow for 
more complicated starting conditions (e.g., holes and filaments; see §4) and multiple disk 
accretion events (e.g., due to merging of subhalos, which supply additional gas; see §5). 
This work makes additional assumptions that can be studied in future work. We have not 
taken into account adiabatic contraction of the dark matter halo. This effect will make the 
gravitational potential well deeper and could be taken into account by changing the values 
of ^0 in the orbit solutions. A deeper potential well leads to greater infall velocities. Since 
the perpendicular component of the velocity is dissipated at the disk surface, the energy 
loss is greater for deeper potentials. For a given angular momentum, a deeper potential 
leads to a smaller centrifugal radius for a given orbit; notice, however, that the angular 
momentum of later orbits may be larger as well (e.g., Elmegreen et al. 2005) and lead to 
a larger centrifugal radius. This current treatment also makes the approximation of zero 
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energy orbits and considers orbital solutions in the inner part of the potential. The latter 
approximation is not severe, since the orbits must be nearly radial in the outer part of the 
potential in order to reach the disk. Similarly, if the initial energy is small but nonzero, the 
orbit will be largely unchanged (where the basic constraint ise=|£^|/^o^l)- 

In addition, the calculation thus far determines the surface density profile resulting from 
the initial infall, i.e., the incoming gas parcels are assumed to stay at the radial locations 
where they fall. Over time, the disk will dissipate energy, transfer angular momentum, and 
spread out. Notice that a moderate amount of disk spreading (angular momentum transfer) 
will make the calculated surface density profiles look more like the exponential reference 
profile over the range of intermediate radii (roughly 2-20 kpc for the parameters used 
here) . Notice also that real disk galaxies depart from an exponential law both on the inside 
and on the outside. The presence of galactic bulges causes the brightness profiles of observed 
disk galaxies to lie above that expected for an exponential disk at small radii, 1-2 kpc (see 
Binney & Merrifield 1998 for a detailed discussion of disk/bulge decomposition). On the 
outside, some disk galaxies show evidence for a cutoff radius Rmax that could be identified 
with Rc in this theoretical treatment; typical values are Rmax ~ 20 — 25 kpc (van der Kruit 
& Searle 1981; Wainscot et al. 1988; Barteldrees & Dettmar 1994; Binney & Merrifield 
1998). 



4. GENERALIZED COLLAPSE SOLUTIONS 

4.1. Pcirtial radial coverage — shells and holes 

This subsection considers configurations for which the initial (pre-coUapse) gaseous 
sphere has holes, i.e., an evacuated central region. Suppose that the central portion of 
the sphere, with radius r < Tq, has no gas. The corresponding surface density profile is that 
of the original sphere (with no central hole) with the contribution from a smaller sphere 
(0 < r < To) removed. If the original sphere has radius tt and mass Mt = M[rT), and 
the hole region has radius ro and mass Mq — M(ro), the combined surface density profile is 
given by 

a{zu) = a{zu; Mt) - a{zu; Mq) , (53) 

where the surface density profiles on the right hand side of the equation are those obtained 
via the spherical collapse calculation of the previous section. 

One example of collapse with a centrally evacuated region is shown in Figure 7. Here we 
consider the intermediate mass profile with p = 2 (M ~ r^) with constant rotation rate, and 
remove successively larger holes from the center of the starting state. The four curves shown 
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Fig. 7. — The effects of removing the central regions from the starting mass profile on the 
resulting surface density profiles. The solid curves show the result of a generalized collapse 
with Rc = 24 kpc and with varying amounts of mass removed from the central regions. The 
upper curve shows the baseline result for a complete gaseous sphere. The other three solid 
curves show the surface density profiles when a central fraction / of the mass of the starting 
sphere is removed, where / = 0.25, 0.50, and 0.75 (from top to bottom). The dashed curve 
shows an exponential surface density profile with the same mass as that of the complete 
spherical starting state and with scale length R^^ — 5 kpc. The scale on the vertical axis is 
arbitrary. 
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in the figure correspond to the full sphere (top curve) and increasing fractions of the mass 
evacuated (25% to 75%). For comparison, the dashed curve shows an exponential surface 
density profile with the same mass as the surface density resulting from the collapse of the 
full sphere. As before, the rotation rate is chosen so that Rc — 24 kpc. 



4.2. Partial angular coverage 



We can also use this formalism to determine the surface density profiles resulting from 
pre-collapse states with partial angular coverage. In general, the starting mass distribution 
can be restricted in both azimuthal angle ip and the polar angle 6. In this problem, the 
starting azimuthal angle determines the azimuthal angle where a given gas parcel will join 
the galactic disk, but it otherwise plays no role (e.g., starting angular momentum depends 
on (r, 6) but not (p). Furthermore, after gaseous material has joined the disk, the differential 
rotation within the disk will eventually make the disk axisymmetric, thereby erasing any 
dependence on (p. Thus, if we consider starting states that are restricted in azimuthal angle 
(f (spheres with longitudinal slices removed), the net result is to reduce the total mass in 
the disk by a constant factor, but otherwise leave the surface density unchanged. We thus 
consider the case of restricting the starting state in polar angle 9. 

Suppose we remove the polar cones from the starting state so that the polar angle is 
restricted to the range < ^ < tt — ^m- The surface density contribution at a given radial 
location in the disk is the same as that prodTiccd by the collapse of the inner portion of 
the sphere with mass M(r = vom), where is given by the intersection of the ray 9 — 9m 
and the locus of points for which material falls to radial position zu, i.e., the locus given by 

rs^/r]^ sm 9o ^ zu . (54) 



Consider the case of an initial mass profile with p — 2 (M ~ r^) and constant rotation 
rate. If the starting state were a complete sphere it would have mass Mt and outer (spherical) 
boundary r^. In this case, the polar cones are removed, where 9m defines the opening angle 
of the cones. At a given radius zu in the disk, the surface density is the same as that for 
mass 

M = Mt ^ ^ ^ . (55) 
Rc sm 9 m 



The surface density profile is thus given by 

Mt 
2'kRc^ 



Mt P dx 2 -,1-1/2 
<^ M = 7^ / — F - 1 > 56) 
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Fig. 8. — Diagram showing the geometry for infall with a polar cone removed. The outer 
sohd circular curves denotes the outer boundary of the collapsing region. The solid ray, 
drawn here at 9m = 7r/4, delineates the evacuated region {9 < 9m)- The dotted curve shows 
the locus of starting points that fall to a particular radius within the disk (taken here to be 
w = 10 kpc). The dashed lines connect the starting locations with the final positions. For a 
given radius in the disk, the surface density is the same as that resulting from the collapse 
of a smaller sphere, one with radius given by the intersection of the ray 9 — 9m with the 
locus of starting points; the dot-dashed circular curve shows this effective boundary for the 
radial location zu — 10 kpc in the disk. 
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where the upper limit of the integration is given by 

. (Rc 1 




} 



(57) 



'M 



Figure 9 shows the resulting surface density profiles for the intermediate starting mass 
profile with M(r) ~ and constant rotation rate. The top curve shows the surface density 
profile for the complete sphere. The next three curves illustrate the effect of removing the 
polar cones from the starting configuration, where the opening angle is taken to be 9m — 
7r/6, 7r/4, and 7r/3. For these three cases, the fraction of the total mass (within the same 
outer boundary) contained in the initial configuration, and hence the resulting disk, is fm — 
V3/2, V2/2, and 1/2, respectively. 



If we restrict the angular extent of a rotating sphere to its limiting form, the result is 
a filament, which provides a reference state for the (pre-coUapse) initial conditions. Notice 
that a filament can be defined in a variety of ways. This treatment starts with a spherical 
initial condition, as considered above, and confines the geometry in both polar angle 6 and 
azimuthal angle ip. In the limiting case, the resulting filament has all of its mass along a fixed 
ray defined by = constant and (/? = constant. Further, this discussion is limited to cases 
where the rotation rate Q = constant, so that the filament has some angular momentum 
about the galactic center but it remains intact as it rotates. These filaments are thus highly 
idealized, but provide a useful benchmark for understanding the nature of galactic infall. 
Keep in mind, however, that these filaments are not the large scale structures (also called 
filaments) that connect galaxies to the cosmic web. 

We begin by considering an isothermal initial configuration, so that M(r) ~ r, which 
will lead to a filament with constant mass per unit length. If we confine the starting sphere 
to a (narrow) range of polar angles 9i < < 62, the resulting surface density can be written 
in the form 



where xi = (sin6'i) and X2 = (sin 62)-^^^. Notice that we can confine the initial state in 
azimuthal angle to a narrow range without affecting result (up to the leading coefficient). In 
the hmit that the range of polar angle is also small, 62 — Oi + 59, the integral in the above 
expression can be evaluated to obtain the result 



4.3. The Limiting Case of Filaments 




(58) 




(59) 
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Fig. 9. — The effects of removing tlie polar caps from tlie starting mass profile on the resulting 
surface density profiles. The solid curves show the result of an axisymmetric collapse with 
Rc — 24 kpc. The upper curve shows the result for a complete gaseous sphere with a p 
— 2 mass profile and constant rotation rate. The next three solid curves show the surface 
density profiles for starting states with polar caps removed, where the polar angle of removed 
material 9m — 7r/6, 7r/4, 7r/3. The dashed curve shows an exponential surface density profile 
with the same mass as that of the complete spherical starting state and with scale length 
i?x = 5 kpc. The scale on the vertical axis is arbitrary. 
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where 6 is the polar angle of the filament (so that 9i ^ 9 and 02 — > ^; 5^ — > 0). Next we 
want to write this expression for the surface density in terms of the effective mass per unit 
length of the filament. The mass contained within the boundary angles 9i and $2 is given by 

M(r) = M,(r/r,) (cos 61 - cos 62) . (60) 

In the hmit of interest, where 92 — 0i — 69<^l, the effective mass per unit length F becomes 

r = {Ms/rs)sme{6e) . (61) 
The surface density can be written in the form 

Tt" 

For comparison, we also consider a filament with mass per unit length F with a polar 
angle 6, where the mass in concentrated along a single ray (centered on the origin). Each 
point along the filament, given by its radial coordinate r, has a particular angular momentum 
given by j = r^QsinO. As a result, conservation of mass implies that each point along the 
filament contributes to the surface density at a particular radius zu within the disk, i.e., 

dm = 27izuadzu = Tdr . (63) 

Using the orbit solution from §2, we can relate the starting radial coordinate r to the resulting 
disk radial coordinate zu, and find 

zu ^ sin e{r^/rl), (64) 



where = TgOJ — r^fi/^/2\&o- Using the above expression to evaluate dr/dvo, the surface 
density (see equation [63]) can be written in the form 

^(^^ = wT^^^k^ ' ^^^^ 
which is in agreement with the result of equation (62). 

As in previous collapse solutions, we can define a centrifugal radius for filaments. For 
the case considered here with constant mass per unit length F, this radius takes the form 

i?c = (MT/Fr,)V^sin^, (66) 

where My is the total mass that falls to the disk (from the filament). With this definition, 
the surface density becomes 

The surface density thus takes a power- law form (here with index q — 3/2) and has a sharp 
truncation radius at Rc- Notice that if we integrate the surface density (eq. [67]) over the 
disk, we obtain the expected disk mass Mr, — Mt (mass conservation is satisfied). 
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5. MULTIPLE COLLAPSE EVENTS 

One defining cliaracteristic of tlie current picture of dark matter lialo formation is tliat 
mergers play an important role — small halos come together to forge larger bound structures. 
In most cases, some portion of the gas in the merging units eventually falls into the center 
of the composite structure, so that galactic disk formation must also involve a series of 
"accretion events" . The goal of this section is to incorporate these multiple accretion episodes 
into our description of gaseous disk formation (numerical simulations automatically include 
this process - see below) . Here we assume that the gas contributed by each merging subhalo 
stacks up inside the gravitational potential of the larger composite structure; after cooling, 
the gas falls inward and thereby adds mass to the disk forming at the center. The galactic 
disk is thus assembled from inside-out, consistent with the standard scenario (e.g., Murali 
et al. 2002, Mailer et al. 2006). Keep in mind that the assimilation of a smaller halo does 
not provide a fully formed disk, but rather provides a new source of gas for the composite 
disk. As a result, each accretion event produces a contribution to the surface density of the 
composite disk. Further, this contribution can be described as a disk surface density profile 
analogous to those calculated above. To summarize, this section develops a description of 
disk formation where a composite disk is assembled through the addition of many separately 
accreted disk layers, where each layer is provided by gas from a merging subhalo.^ 

Numerical simulations of hierarchical structure formation incorporate the idea of mul- 
tiple accretion events in a natural way. The addition of mass into a developing halo is often 
described in terms of the halo mass assembly history (MAH), which has been well studied 
(e.g., Avila- Reese et al. 1998; Firmani & Avil- Reese 2000; Wechsler et al. 2002; van den 
Bosch 2002). This body of previous work indicates that the MAH is generally not domi- 
nated by major mergers; instead the MAH consists of a continuing series of smaller-scale 
events that can be described as a smooth aggregation of mass (Murali et al. 2002; Mailer et 
al. 2006). In the (idealized) model developed herein, we consider each of these small-scale 
accretion events to provide a supply of gas that (after infall) results in a constituent surface 
density profile that contributes to the overall composite disk profile. The smooth nature of 
the MAH suggests that we must consider a large number of accretion episodes (each with a 
relatively small baryonic mass). 

The distribution of angular momentum of the gas added during the accretion events 



^Although the net effect of merging halos is to add disk surface density profiles together, we stress that 

disks themselves are not merging. Instead, the merging halos provide a new gas supply for continued infall 
onto the disk. Notice also that we are making an assumption of timing, i.e., that incoming halos have not 
yet formed their disks. 
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determines the distribution of centrifugal radii for the constituent disks is thus an important 
quantity. We expect the newly added gas to have spin angular momentum inherited from the 
rotation of the subhalo that supplied the gas. In addition, the subhalo and its gas can have 
orbital angular momentum. Numerical simulations indicate that galactic spin parameters 
Xj — and hence the spin angular momentum of subhalos and their g drawn from 

a log-normal distribution (e.g., Bullock et al. 2001); these numerical studies also suggest 
that the width of the distribution a ~ 0.5 — 0.6. Although the angular momentum will 
include an additional orbital component, the latter cannot dominate; otherwise, the total 
angular momentum of the disks would be too large (see below) . The distribution of orbital 
angular momentum is not known, but it will tend toward a log-normal distribution if a 
large number of physical variables play a role in setting the orbits (due to the central limit 
theorem; Richtmyer 1978). Another issue is that the angular momentum vectors (both spin 
and orbital) of the constituents are not necessarily aligned. Since the angular momenta 
add as vectors, this complication will act to make the total (composite) angular momentum 
smaller (than the result of adding scalars); further, this process adds another independent 
physical variable which (through the central limit theorem) will help simplify the composite 
distribution. As a result, we expect the distribution of total angular momentum to have 
nearly a log-normal form. 

Before proceeding further, it is useful to summarize the ingredients of the calculation. 
The vahdity of this approach relies on the following assumptions: [1] The composite galactic 
disk forms as as superposition of gaseous layers, where each layer can be considered as a 
constituent disk with a power-law surface density distribution. We first consider equal mass 
constituents and then generalize to varying disk masses (which can be either correlated or 
uncorrelated with disk radius). [2] Each gaseous layer is assembled from the gas contributed 
by a subhalo that merges with the growing galactic halo. [3] During the merging process, 
the gas of each "accretion event" retains its component of angular momentum aligned with 
the overall sense of rotation of the galaxy. [4] During the course of a particular accretion 
event, inf ailing gas conserves its angular momentum and produces a surface density profile 
in accordance with the solutions of §2 - 4. The centrifugal radius of each disk layer is then 
determined by the total (spin and orbital) angular momentum of its progenitor subhalo. 
[5] Over the course of multiple merger events, the resulting distribution of (total) angular 
momentum is log-normal; as a result, the distribution of centrifugal radii is also log-normal 
with a well-defined peak value Rq and width a. 

This calculation begins with the simplest case of equal mass constituent disks with the 
same general form, but with different centrifugal radii (which correspond to to different 
angular momenta of the pre-coUapse gas). The basic problem is thus to add together an 
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ensemble of contributing disk profiles with the form 




(68) 



where x = w / Rj and Rj is the centrifugal radius of the jth constituent disk (from the 
jth accretion or merger event). The function g{x) defines the shape of the edge, so that 
g{l) 0. The analysis given above shows that typical edge functions are g{x) = cos^^ x 
and g{x) = {1 — x^Y^^ . The integer N represents the number of events, and hence the total 
mass of the composite disk is specified by Mq times a dimensionless factor given by the 
integral / = 27r x^~'^g{x)dx. 

In our formalism, the above considerations imply that the ensemble of centrifugal radii 
are drawn from a log-normal distribution. Let Ro be the radius at the peak of the distribution 
and let a denote the width of the distribution. The composite disk surface density then takes 
the form 



The sum can be evaluated by random sampling to determine the composite surface density 
profile of the disk. One example is shown in Figure 10, which uses the surface density 
distribution of equation (28) to specify the constituent profiles. The sum (69) is evaluated 
in the large N limit. As expected, the composite surface density profile depends on the 
width of the log-normal distribution. For small values of the distribution width a, all of 
the components are nearly the same and the composite profile looks much like that of the 
original constituents. In the limit of large width a, the composite profile approaches a power- 
law form (fj ~ zu~^, as derived below) over the range of radii of interest. In the realm of 
intermediate a, the composite surface density profile is close to the exponential reference 
profile (shown by the dashed curve in Fig. 10). All of these results are robust (in that 
they are largely independent of the assumed constituent profiles) and can be understood 
analytically, as discussed below. 

To study the problem of constructing composite surface density profiles, we need to 
evaluate the sum in equation (69). We consider the limiting case in which the disk edges are 
"sharp", so that the edge function is a step function g{x) — Q{l — x). For a given reafization 

of the disk parameters, we then order the disk radii Rj, from largest to smallest, and define 
F{zu) to be the fraction of constituent disks with w < Rj. The sum then becomes 




(69) 




(70) 
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Fig. 10. — The composite disk surface density profile resulting from the addition of a large 
number of constituent disks. Each component disk has the form a ~ cos~^(ro/i?j), i.e., a 
power-law with a well-defined outer edge at the centrifugal radius Ry The radii Rj are drawn 
from a log- normal distribution with varying dimensionless widths = 0.1 (lowest curve), 
0.5, 1.0, and 1.5 (top curve). The radial scale is set by the peak of the distribution at i?o- 
The dashed curve shows a reference exponential profile, where the scale length R^ = 0.35i?o- 
Notice that the intermediate value of o" = 0.5 comes close to producing an exponential profile. 
Notice also that the scale on the vertical axis is arbitrary. 
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Since the Rj are log-normally distributed, we can convert the sum to an integral in the limit 
of large N, where a large number of constituent disks are being assembled: 

aT{w)=MoW-m,'^'-''> r ^ (71) 

Jzo{-cj) v27r(T 

where zq — ln{zu/RQ) and z — \n{Rj/Ro) is the natural logarithm of the centrifugal radii 
being integrated over. Recall that Ro is a radial scale corresponding to the peak of the 
log- normal distribution and a is the distribution width. Because the sum (eq. [70]) only 
includes disks with outer radii Rj > zu, the corresponding integral (eq. [71]) has a limited 
range of integration from Zq{zu) out to z ^ oo. This integral can be written in terms of error 
integrals, which have a number of known properties (Abramowitz & Stegun 1972; hereafter 
AS72). In particular, we find 

ariru) = Motz7-«i?-^'-'^ g^^-^^'^'/a^ / ^^e"^' , (72) 

where = [-^o + (2 — q)a'^]/{V^a). To approximate the integral, we make use of the well- 
known inequality, 

1 f°° 1 

< e«' / e-''dt < / , (73) 



which is valid for ^ > (AS72). This constraint allows us to write the composite surface 
density profile in the form 



aTizu) ^ , , (74) 

where both zq and are functions of the disk radial coordinate vo. The constant Ce must 

lie in the range 4/7r < Cc < 2, and can be chosen to minimize the error. If we take Cg 
= 3/2, the approximation to the integral in equation (73) is accurate to ~ 10% for all 
allowed values of the parameters, and is accurate to ~ 2% over the range of values needed 
to evaluate the surface density profile. Notice that this form is only valid for > 0, which 
implies w > i?oexp[— (2 — q)a^]. For sufficiently large a, however, this form is valid over 
most of the radial range of interest. These trends are illustrated in Figure 11, which shows 
the composite surface density calculated from random sampling compared to that obtained 
using the approximation of equation (74). 

Now consider the composite distribution in the limiting cases. In the limit where the 
width a of the distribution is narrow, the normal distribution (in the logarithm) becomes a 
Dirac delta function (see eq. [71]) and we recover the expected result 

--^(^)'«(^-/^")- (^«) 
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Fig. 11. — The composite disk surface density distribution obtained from random sampling 
compared to that obtained from the asymptotic approximation derived in the text. The 
sohd curve shows the composite profile obtained by sampling = 10^ constituent surface 
density profiles, with simple power-law form a ~ l/c^?, truncated at centrifugal radii Rj, 
which arc sampled from a log-normal distribution of width a = 0.5 (the radial scale is set by 
the peak Rq of the distribution of centrifugal radii). The dotted curve shows the composite 
surface density obtained from the approximation of equation (74), which provides a good 
representation for r ^ 0.6i?o- The dashed curve shows an exponential profile for comparison. 
The scale on the vertical axis is arbitrary. 
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In the opposite limit of large width o", equation (74) is valid for essentially all radii as long 
as < 2, and we find that the composite surface density takes a power-law form, i.e., 



Mn 



27rCT(2 - q) 



-VJ 



(76) 



This resulting power-law form aT oc w~'^ is independent of the index q of the constituent 
disk profiles, as long as g < 2. Notice that the results of §3 show that g < 2 for essentially all 
cases of physical interest, as indicated by equation (44). For large values of the initial index 
0' > 2, all of the mass is concentrated toward the galactic center and the composite surface 
density profile always takes the form ut oc This result follows from equation (72). As 
a — > GO, ^0 becomes large and negative, the integral in equation (72) approaches a constant 
value (-\/7r)) and the composite surface density takes the form ut ~ vj~'^. 

In the limiting case of o" — > 0, the disk edge becomes sharp, which (in a sense) is much 
steeper than an exponential falloff (at large vo). In the opposite limit of large o", the disk has a 
power-law form ~ va~'^, which falls off more slowly than an exponential function. One might 
suspect that for some intermediate value of a, the disk would obtain a roughly exponential 

shape, at least for some range of radii. This behavior docs in fact occur, as shown in Figure 
10, and can be understood in terms of the analytic approximation derived above. For radii 
sufficiently far from the origin so that equation (74) is valid, we can define an effective scale 
length ^cr for the surface density profile via 



1 dai 



(Tt dzu ZD 



1 I ^ ^0 

2 + -^ + 



(^o + (2-?)a2)' + 2a2ce] 



-1/2 



(77) 



If this scale length laizu) were exactly constant, when considered as a function of ro, the 
surface density profile would be exponential. In practice, this function is slowly varying. 
Specifically, over the range of radii \ < w jR^ <^ (which corresponds to < ^ 1.609), 
the scale length 1^ varies by 10-12% for q in the range 0.5 - 1.5, where the width a is chosen 
to minimize the variation. This optimization requires values of cr = 0.3 - 0.5, which are 
comparable to, but somewhat smaller than, the widths of the distributions deduced from 
numerical simulations (where a fa 0.5 — 0.6; Bullock et al. 2001). 

The discussion thus far assumes that each constituent disk has a constant mass (equal to 
Mq/N). However, the formalism is sufficiently robust to consider the general case of varying 
disk masses. The composite profile will not change as long as the disk masses are uncorrelated 
with the centrifugal radius Rj of the individual components. Specifically, suppose that the 
mass of each component disk has the form 



M^^Mofj{Rj/Ro)\ 



(78) 
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where fj is a random variable with unit mean. The masses are uncorrelated when a — 
and show a correlation with the centrifugal radius for a 7^ 0. If we propagate this ansatz 
through the formalism given above, the composite disk profile takes the form 

Mo e~^o/2^' 

where the variable .^0 = [-^o + (2 — — q;)(T^]/(^/2o"). Notice that the effect of the random 
variable fj is completely washed out in the limit of a large number of accretion events. In the 
limit of a narrow distribution, we again recover the expected result that the composite distri- 
bution is the same as that of the components (eq. [75]). In the limit of a wide distribution, 
the composite surface density profile becomes a power- law with index 2 — a, independent of 
the index q of the component profiles, as long as q < 2 — a. For "large" values oi q > 2 — a, 
the composite surface density takes the form of the original profiles so that ar ~ W'^. 

The results presented thus far have been calculated in the limit of large N, so that the 
surface density of the composite system approaches a particular smooth form. An interesting 
question is to ask how large N must be in order for the problem to be safely in the large N 
limit. One way to answer this question is to calculate the RMS deviation of the difference 
between the composite surface density obtained with a "small" value of and that of the 
large ^ 00 limit. Notice, however, that for small A^, each realization of the random values 
(the Rj) will lead to a somewhat different profile and hence a different measure of the RMS 
error. Roughly speaking, for A^ = 10, we get deviations of about 50 percent, for A^ = 100 the 
deviations are about 10 percent, and for A^ = 1000 the deviations are only about 2 percent. 
These deviations are computed over the range < uj < 2Rq] they are greater at larger radii. 
Another way to address this issue is to plot the resulting surface density profiles for varying 
numbers A^ of the constituent disks. The result is shown in Figure 12 for logarithmically 
spaced A^ = 1, 4, 16, 64, and 256. As the number of constituent disks A^ increases, the 
profiles approach that of the asymptotic (large A^) limit shown by the heavy solid curve in 
Figure 12. Keep in mind that the profiles will vary greatly from realization to realization for 
small values of A^; in other words, the small- A^ profiles shown in the figure will look different 
every time they are computed (with independent sampling). Nonetheless, for A'> 16, the 
composite surface density has roughly the same form as that of the asymptotic limit out to 
r ^2Rq. Similarly, for A^> 64, the surface density profile will be similar to the asymptotic 
hmit out to r ~ 3i?o- This result is important, since galaxy sized halos are thought to be 
assembled from moderate numbers of component halos (Wechsler et al. 2002). 

This formalism can be applied to any scenario in which the net effect is to add surface 
density profiles together. As one example, it can account for cosmological evolution of the 
background potential. In the approximation used here, the surface density profiles of the 
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Fig. 12. — The composite disk surface density profile for varying number N of constituent 
disks. Each component disk has the form a ~ zu~^ cos~^ {lu / Rj) , where the radii Rj are 
drawn from a log-normal distribution with dimensionless width a = 0.5 (and peak at Rq). 
The heavy solid curve shows the resulting profile in the limit of large A^; the lighter solid 
curves show the results for = 1, 4, 16, 64, and 256, where the curves approach that of 
the large limit with increasing A^. Notice that for small A^, the profile will vary from 
realization to realization due to incomplete sampling of the distribution of disk radii. The 
dashed curve shows an exponential surface density profile as a reference point. The scale on 
the vertical axis is arbitrary. 
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constituent disks add together so that the final state does not depend on the order in which 
they are assembled. In a cosmological setting, however, dark matter halos become steadily 
larger, both from continued infall from large distances and through the mergers of smaller 
halos. As a result, the depth of the gravitational potential well (\E'o) increases with time. 
If the constituent disks result from the collapse of gas that attains hydrostatic equilibrium 
(§3.5), then the pre-coUapse extent of the gas will be smaller at later times (for a fixed mass), 
and the corresponding centrifugal radius Rc will be a decreasing function of time if angular 
momentum is constant. On the other hand, some observations (Elmegreen et al. 2005) 
and theoretical work (Firmani & Avila-Reese 2000) indicate that disk radii should increase 
with time. This result can be made consistent with our multiple accretion scenario if the 
specific angular momentum increases for subsequent accretion events. In any case, the time 
dependence of the centrifugal radii can be taken into account in the distribution of Rj. If the 
distribution is log-normal, as assumed here, then the final result depends only on the peak 
value Ro and width & of the distribution. However, this procedure can be generalized to 
other distributions. As another example, some numerical simulations suggest that baryons 
can transfer angular momentum to the dark matter so that gas falling in later would have a 
systematic shift in its angular momentum distribution (private communication, anonymous 
referee) ; the effect of this process on the resulting surface density profile of the disk could be 
studied with the formulation developed here. 

This analysis shows that a composite disk can display a nearly exponential surface 
density profile, provided that the distribution of centrifugal radii has an intermediate width 
{a 0.5) and the constituent disks have power-law indices q < 2. Both of these conditions 
are expected to be met, as outlined above. Notice that when an exponential surface density 
profile is realized, the scale length is roughly R^ ~ 0.35i?o (Figs. 10, 11, 12), where Ro is 
the centrifugal radius at the peak of the distribution. To produce scales i?x ~ 4 — 5 kpc, as 
indicated by observations (Binney & Merrificld 1998), the radius Ro must lie in the range 
i?o = 11 ~ 14 kpc. The centrifugal radius of an individual constituent disk depends on the 
spin parameter A through a relation of the form Rc — aXrg, where et is a dimensionless 
parameter of order unity. For example, if the halo has a Hernquist profile, the gas density 
follows that of the dark matter, and all of the gas out to r = falls in, then a ^ A. With 
typical halo scale radii = 50 - 70 kpc and spin parameters A = 0.035 - 0.04 (at the peak 
of the distribution - see Bullock et al. 2001 and Wechsler ct al. 2002), the expected peak 
value of the centrifugal radius would lie in the range i?o = 7 - 11 kpc. These values of A are 
thus comparable to, but somewhat lower than, those required to produce the correct disk 
scale length. However, the values quoted above only include the spin angular momentum 
of the constituent gas modules. As these subunits fall into the larger halo and merge, and 
thereby add their gas to the growing structure, they contain orbital angular momentum 
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in addition to that described by the spin parameter A. If the orbital angular momentum 
is roughly comparable to (but smaller than) the spin angular momentum, the constituent 
disks will have the proper total angular momentum for this multiple infall process to produce 
exponential composite disks with the observed scale length. 

6. CONCLUSION 

6.1. Summciry of results 

This paper has studied the assembly of galactic gaseous disks by following orbit solutions 
for baryonic gas parcels as they fall through the potential wells of dark matter halos. Our 
general analytic formulation of the problem (as presented in §2) is robust and can be used 
in variety of contexts. In this paper, we calculate the velocity fields of the incoming gas, 
the density distribution of the gaseous material as it falls inward, and the surface density 
profiles of the resulting gaseous disks. Our specific results are summarized below: 

[1] The generic outcome for a gaseous disk produced by a single spherical collapse is a 
power-law surface density that is truncated at a centrifugal radius Rc (see Figs. 2, 4, and 
5). Here, the centrifugal radius Rc is the radius to where the material with the highest 
(initial) specific angular momentum falls during the collapse. This result suggests that 
the "fundamental units" of galactic disk assembly arc not exponential disks, but instead 
are truncated power-laws. For a constant rotation rate of the initial gaseous sphere, the 
power-law index q of the surface density distribution {a oc zu^'^) depends on the initial mass 
distribution such that q = 2 — p, where p is the power-law index of the starting mass profile 
(M oc r^). For a rotation rate that varies with radius according to ^2 ~ r^*, the power-law 
index of the resulting disk surface density is given byg = 2 — p/(2 — s). We also determined 
the density profiles appropriate for hydrostatic equilibrium of both isothermal and polytropic 
gas in the background potential of a dark matter halo (§3.5) and followed the subsequent 
collapse (Figs. 5 and 6). For polytropic equations of state with relatively large indices n, 
the resulting surface density profiles are close to exponential (except near the centrifugal 
barrier) . 

[2] In addition to determining the disk surface density profiles, this calculation provides 

an analytic description of the dynamics that builds up galactic disks. This description 
includes the velocity fields of the gas during the collapse (eq. [8]), the density distribution 
of the infalling baryonic material (eqs. [10] and [16]), and the mapping between the starting 
locations in the prc-collapsc cloud and the "final" radial positions within the disk (e.g., Figs. 
3 and 8). We also obtain analytic forms for the hydrostatic equilibrium profiles of polytropic 
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gaseous spheres in the gravitational potential of a dark matter halo (eqs. [45] and [49]). 

[3] This collapse picture has been generalized to include less idealized starting configura- 
tions (§4). The first generahzation considers the starting gaseous spheres to have evacuated 
regions in their centers. The signature of these central holes in the starting states is to 
produce impoverished central regions in the corresponding gaseous disks (see Fig. 7). The 
starting states can also be considered to have limited angular extent. The removal of the 
polar caps from the initial states also leads to a deficit of material in the disk centers (relative 
to the spherical initial state) as shown in Figure 9. 

[4] This formulation of the collapse problem can be further generalized to include multi- 
ple accretion events, which are expected to occur due to merging of individual halos to form 
larger, composite structures (§5). Each incoming halo provides a new supply of gas, which 
subsequently cools and falls inward to produce another component to the disk. Here we 
consider the addition of a large number of such constituent disks, where the centrifugal radii 
are drawn from a log-normal distribution (consistent with the numerical finding that halo 
spin parameters display a log-normal distribution) and find an interesting and somewhat un- 
expected result: The addition of a large number of power-law disk profiles with index q <2 
results in a composite surface density profile of the form a ~ vd~'^ in the limit of a wide 
distribution for Rj, where this result is independent of the index q (provided that ^ < 2). 
In the limit of a narrow distribution of the composite surface density is the same as 
the constituent surface density profiles. In the intermediate regime, however, the composite 
surface density can display a nearly exponential form over a range of radii (Figs. 10 and 11). 
The asymptotic forms are obtained in the limit of large numbers N of constituent disks, and 
the approach to large N is shown in Figure 12. Although the large N hmit is approached 
rather slowly, this sequence shows that tens of subunits (constituent disks) are necessary to 
produce a smooth disk with a nearly exponential surface density profile. 



6.2. The quest for exponential disks 

One of the defining features of observed disk galaxies is that their surface brightness 
distributions display a nearly exponential form. One might hope that numerical simula- 
tions of galaxy formation would naturally produce gaseous disks with exponential surface 
density profiles, so that star formation could take place in situ and thereby account for the 
observations. Unfortunately, the production of exponential surface density profiles is not au- 
tomatic. Building upon earlier work (Dalcanton et al. 1997, van den Bosch 2001, Kaufmann 
et al. 2006, and many others), this paper shows that the fundamental disks produced by 
a single collapse event have the form of a truncated power-law rather than an exponential 
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profile (Figs. 2 and 4). The transformation between these truncated power-law forms and 
exponential profiles most likely involves several processes, as outlined below: 

[A] The observed disk brightness profiles are close to exponential only over a hmited 
range of radii (roughly 2-20 kpc) and a correspondingly small range of scale lengths (roughly 
r/i?x ~ 0.5 — 5). The inner regions depart from an exponential form due to contamination 
from bulge components and the outer regions are too faint to be observed. Further, as 
shown by Figure 2, the difference between a truncated power-law and an exponential is not 
too large over a limited range of radii. One of the largest differences between the power- 
law profiles and an exponential profile occurs at small radii, where the power-law profiles 
have cusps; this discrepancy is largely alleviated, or at least masked, by the presence of the 
bulge component (for r < 2 kpc). In addition, if the initial (pre-coUapse) gas departs from 
a complete spherical form, either through inner evacuated regions or the removal of polar 
caps, then the resulting surface density profiles will not have inner cusps (Figs. 7 and 9), or 
they will be less pronounced. 

[B] The calculation presented in this paper accounts for only the formation of the disk, 
not its subsequent evolution. If galactic disks are produced with sharp edges (as indicated 
here for the case of a single collapse event), then gravitational instabilities will grow robustly 
and cause the disk to spread. Although previous authors (Dalcanton et al. 1997) have 
argued that the time scale for such instability is too long, the quoted (long) time scales are 
appropriate for disks that already have soft (e.g., exponential) edges. If disks have sufficiently 
hard edges, as produced here, then gravitational instabilities can grow on the dynamical time 
scale of the outer disk edge (Adams et al. 1989; Ostriker et al. 1992). This time scale is given 
by the inverse of the orbital frequency at the centrifugal radius, roughly r = Rc/vrot ~ 100 
Myr {Rc/20 kpc). This time scale is still somewhat long to account for all of the necessary 
angular momentum redistribution, but gravitational instabilities will have time to grow and 
act to spread out the disk. Notice that a truncated power-law profile that spreads out will 
become more like an exponential profile. 

[C] Exponential profiles are observed for the surface brightness distributions of galactic 
disks, not the surface density of the gas, so it is possible for the star formation process to 
help produce the observed exponential forms. A good demonstration of this effect is found 
in Lin & Pringle (1987), where initial surface density profiles roughly like those of Figure 
2 spread out under the action of disk viscosity and convert some fraction of the mass into 
stars. The resulting surface density of the stars is much more like an exponential form than 
the starting gas profile. 

[D] Finally, as shown in §5, the addition of multiple disk surface density profiles can 
produce a composite profile that displays a nearly exponential form — much more like an 
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exponential than the initial (constituent) starting profiles. The ingredients necessary to pro- 
duce a nearly exponential profile are a large number of constituent disks and a log-normal 
distribution for their centrifugal radii. This result is independent of the shape of the con- 
stituent disk profiles as long as they have effective power-law indices q < 2. Since the current 
paradigm of structure formation indicates that galaxies are assembled through a process of 
hierarchical merging, multiple accretion events are highly likely. Further, numerical simula- 
tions suggest that the distribution of halo angular momentum is log-normal (Bullock et al. 
2001) so that the distribution of centrifugal radii is also expected to be log- normal. 

In practice, all four of these considerations (and perhaps others) will conspire to produce 
the observed exponential disk profiles. Individual infall events will each produce disk surface 
density profiles that are truncated power-laws, but multiple infall events will jointly build 
up a nearly exponential disk (as shown in Fig. 10). The effects of star formation, viscous 
spreading, and disk evolution via gravitational torques will act to smooth out the disk 
and drive the surface density towards a more exponential form. Given the presence of the 
bulge contribution at small radii and the limited dynamic range of the observed profiles, the 
exponential disks in observed galaxies can be understood through the joint action of these 
processes. 
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